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Are gauge shocks really shocks? 
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The existence of gauge pathologies associated with the Bona-Masso family of generalized harmonic 
slicing conditions is proven for the case of simple 1+1 relativity. It is shown that these gauge 
pathologies are true shocks in the sense that the characteristic lines associated with the propagation 
of the gauge cross, which implies that the name “gauge shock” usually given to such pathologies is 
indeed correct. These gauge shocks are associated with places where the spatial hypersurfaces that 
determine the foliation of spacetime become non-smooth. 
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I. INTRODUCTION 

When studying the dynamical evolution of spacetime, 
it is important to choose coordinates that allow one to 
cover as large a region of spacetime as possible with¬ 
out becoming pathological. In the 3+1 approach, this 
choice of coordinates naturally separates in two differ¬ 
ent aspects: The choice of a time coordinate, that is, of 
a foliation of spacetime into spatial hyper-surfaces (also 
known as the “slicing”), associated with the lapse func¬ 
tion a , and the choice of the way in which the spatial 
coordinates (the “time lines”) propagate from one hyper¬ 
surface to the next, associated with the shift vector (3 l . 

With respect to the choice of slicing, many different 
ways to choose the lapse function are possible. One can, 
for example, specify the lapse directly as a function of the 
geometric variables. Alternatively, one can obtain the 
lapse as a solution of some specific differential equation. 
Elliptic slicing conditions are typically obtained when one 
enforces some geometric condition on the spatial hyper¬ 
surfaces. An example of an elliptic slicing condition is 
“maximal slicing” [jj, which requires that the spatial vol¬ 
ume elements remain constant during the evolution and 
has strong singularity avoiding properties, making it par¬ 
ticularly well suited for studies of black hole spacetimes. 

Another possibility is to specify an evolution equation 
for the lapse function and evolve this function in time 
as just another dynamical quantity. This last approach 
has the advantage that it is much easier computation¬ 
ally to evolve the lapse than to solve an elliptic equation 
for it. One particular family of evolution type slicing 
conditions was introduced by Bona and Masso as a gen¬ 
eralization of the harmonic time coordinate condition |2j. 
This family has the important property of allowing one 
to construct strongly hyperbolic formulations of the Ein¬ 
stein evolution equations that include the slicing condi¬ 
tion. Also, some members of the Bona-Masso family have 
been shown to mimic the singularity avoiding properties 
of maximal slicing. 

In 1997,1 studied different members of the Bona-Masso 
family and found that, unless a specific condition was im¬ 
posed, gauge pathologies could easily develop Q. I called 
those pathologies “gauge shocks” as they appeared as 
discontinuities in the lapse and metric functions that de¬ 


veloped from smooth initial data. In Q I strengthened 
the case for the existence of these gauge shocks by pro¬ 
viding a purely kinematic argument independent of the 
Einstein equations, and more recently in |5} gauge shocks 
have again been predicted together with blow-ups asso¬ 
ciated with constraint violation. Still, there have been 
some doubts in the numerical relativity community about 
both the relevance and the reality of gauge shocks. Re¬ 
cently it has even been claimed by Bona et al. Tij that 
genuine gauge shocks are completely discarded. 

The purpose of this paper is therefore to provide a 
rigorous proof in the case of simple 1+1 dimensional rel¬ 
ativity that gauge pathologies associated with the Bona- 
Masso family of slicing conditions do develop, and that 
these pathologies are genuine shocks in the following 
sense: They correspond to discontinuities in the solu¬ 
tions of the underlying hyperbolic differential equations 
that develop from smooth initial data and are associated 
with the crossing of the characteristic lines. 

It is important to stress the fact that the term “shock” 
will be used here strictly in this restricted sense of the 
crossing of the characteristic lines. In hydrodynamics, 
the word shock is associated also with the subsequent 
propagation of these discontinuities (“shock waves”). 
However, since at a discontinuity the differential equa¬ 
tions break down, one can only talk about so-called 
“weak” solutions. Such weak solutions are not unique, 
and one needs to specify an extra physical principle 
known as an “entropy condition” to identify the correct 
solution. In the case of gauge shocks, however, once the 
discontinuity develops our gauge has in fact broken down. 
It is unlikely that one can find an entropy condition in 
this case as the gauge can be chosen arbitrarily, so any 
weak solution should be acceptable in principle (though 
they would all be singular). Notice also that we are re¬ 
ally interested in finding how such gauge shocks can best 
be avoided since discontinuities in the gauge variables are 
clearly undesirable. We are not interested in propagating 
these discontinuities. 0 

This paper is organized as follows. In Sec. Eli review 
the basic equations of 1+1 relativity, and show that they 
can in fact be written as a system of conservation laws, 
both in terms of the geometric variables and in terms of 
the eigenfields. My main argument for the existence of 
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gauge shocks is presented in Sec. m first showing that 
blow ups of the eigenfields can develop in finite time, and 
then showing that the characteristic lines associated to 
the gauge propagation do cross. I also show a numerical 
example showing how a guage shock develops. I conclude 
in Sec ITvl Appendix El describes the initial data used 
for the numerical simulations, and finally Appendix m 
counters the argument given by Bona et al. in :6j to 
claim that gauge shocks do not occur. 


II. EINSTEIN EQUATIONS IN 1+1 
DIMENSIONS 

Let us consider vacuum general relativity in one spatial 
dimension. It is well known that in such a case the grav¬ 
itational held is trivial and there are no true dynamics. 
However, one can still have nontrivial gauge dynamics 
that can be used as a simple example of the type of be¬ 
havior one can expect in the higher dimensional case. 

As slicing condition we will use the Bona-Masso family 
of generalized harmonic slicing conditions j2| 

d t a = -a 2 f(a)I< , (2.1) 

with K = K* the trace of the extrinsic curvature. 

The ADM J3, .8] evolution equations in the 1+1 case 
can be written in first order form as 

dta = - a 2 fl <, ( 2 . 2 ) 

d t g = -2 agK, (2.3) 

and 

d t D a +d x (afK) = 0, (2.4) 

d t D g + d x {2aK) = 0, (2.5) 

d t K + 8 X ( aDJg) = a (A 2 - D a D g /2g) , (2.6) 

where we have defined g g xx , D a := d x In a and 
D g := d x In g. 

Before doing an analysis of the characteristic structure 
of this system of equations, it is important to notice that 
the evolution equation for A can in fact be rewritten as 
a conservation law in the following way 

d t (> /2 A) + d x ( aD a /g V 2 ) = 0 . (2.7) 

If we now define the vector v := (D a , D g , A), with 
K := g 1 / 2 A, then the evolution equations for the first 
order variables can be written as a conservative system 
of the form 


d t v + d x (M v) = 0 , (2.8) 

with the characteristic matrix M given by: 

/ 0 0 af/g 1 / 2 \ 

M = I 0 0 2a/g 1 / 2 . (2.9) 

\ «/ff 1/2 0 0 / 


The characteristic matrix has the following eigenvalues 


Ao = 0 , A± = ±a(f/g) 1/2 , 

(2.10) 

with corresponding eigenvectors 


3) = (0,1,0), e±= (/,2,+Z 1 / 2 ) . 

(2.11) 

Since the eigenvalues are real for / > 0 and the eigenvec¬ 
tors are linearly independent, the system 12.811 is strongly 
hyperbolic. The eigenfunctions are given by 

LU = R l v , 

(2.12) 

with R the matrix of column eigenvectors. We find (using 
an adequate choice of normalization, see below) 

Wo = D a /f - D g / 2 , w± = A ± A*// 1 / 2 , 

(2.13) 

which can be easily inverted to give 


~ (w++u;_) 

K = 2 ’ 

(2.14) 

= I 1 ' 2 (u + -u-) 
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(2.15) 

_ (w+-w_) 0 _. 

D g - /2 2 uj 0 ■ 

(2.16) 


It is important to notice that with the eigenfunctions 
scaled as above, their evolution equations also turn out 
to be conservative and have the simple form: 


d t ti + d x (A w) = 0 , (2.17) 

with A = diagjAi}. If, however, the eigenfunctions are 
rescaled in the way w' = Fi(a,g) u>i , then the evolution 
equations for the will in general no longer be con¬ 
servative and non-trivial sources will be present. The 
important point is that there is in fact one normalization 
in which the equations are conservative, namely the one 
given in 12.1311 . 


III. GAUGE SHOCKS IN 1+1 RELATIVITY 


There are two different ways in which one can see that 
the evolution equations derived in the previous section 
develop singular solutions. Let us start by looking at the 
evolution equations for the traveling eigenfunctions u>± 

d t u)± + d x (A±w±) = 0 . (3.1) 


We now rewrite these equations as 

dt u >± + A ±d x uj± = —u±d x \± . 


(3.2) 


Using the expressions for A±, and denoting f = df /da, 
one finds 


_ t «/ 1/2 a +/ 1/2 
9 x X± — +— , /n d x g ± 


2<? 3 / 2 


9 


1/2 


1 + 


Of 

2 / 


d x a 


= A± 
= A± 


of A Da _ Dg_ 
1 2 ) f 2 

1 2 ) 2/!/2 


- U! 0 


, (3.3) 
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and finally 


dtu>± + \±d x io± — 


A±w± 



af'\ lo + — 

~2~ J 2/V2 W ° 


(3-4) 


Assume now that we are in a region such that 
too = LO- =0. It is clear that ujq will not be excited, since 
it does not evolve, nor will lo_ be excited, since from the 
equation above we see all that its sources vanish. The 
evolution equation for a>+ then simplifies to 

d t w + + A +d x w+ = ^/2 (l - / - . (3.5) 

This equation shows that w .(_ will blow up along its char¬ 
acteristics unless the term in parenthesis vanishes: 


Eigenfield co + 



(X f' 

1-/-^T = 0. (3.6) 

The last condition has been derived several times be¬ 
fore III and can be easily integrated to give 

f(a) = 1 + k/a 2 , (3.7) 


with k an arbitrary constant. Notice that harmonic slic¬ 
ing given by /= 1 is of this form, but / = constant ^ 1 is 
not. Notice also that “1+log” slicing for which / = 2/a, 
though not of the form in . nevertheless satisfies con¬ 
dition (EDI at places where a = 1. 

In some cases it is even possible to predict exactly when 
a blow up will occur. In order to see this we first define 
the rescaled eigenfunctions := au>±/g 1 ^ 2 . For their 
evolution equations we now find 


dt^l± + Ai^ri-i- — 


+ 



an n| 

2/2 


af \ n±n T 
1 2 ) 2 


(3.8) 


Notice that with this new scaling, all contributions from 
loq to the sources have disappeared. If we now assume 
that we have initial data such that = 0, then the 
evolution equation for f l + simplifies to 




/ 


af \ 

2/2 


(3.9) 


which can be rewritten as 

dV+ = (,_ f a f\ 

dt V 1 2/2 


(3.10) 


with d/dt the derivative along the characteristic. It is 
clear that we have the same condition for avoiding blow¬ 
ups as before. But the last equation has a very important 
property: for constant / the coefficient of the quadratic 


FIG. 1: Evolution of the eigenfunction u>+ for initial data rep¬ 
resenting a pulse traveling to the right. The dotted line shows 
the initial data, and the solid line the numerical solution at 
t = 10. Notice how u>+ has developed a large spike. 


source term is itself also constant. In that case the equa¬ 
tion can be easily integrated to find (assuming f ^ 1) 


= 




i - (i - /) t/2 


(3.11) 


where = Q + (t = 0). The solution will then blow up 
at a finite time given by t* = 2/[(l — /) S2(|_]. Clearly, this 
time will be in the future if (1 — /) Q+ > 0, otherwise 
it will be in the past. Since in general will not be 
constant in space, the first blow-up will occur at the time 


T* = 2/[(l - /) max(fi° (*))] . (3.12) 

Figures 0 and Q show the numerical evolution of the 
eigenfield u>+ and the lapse function a, in a case when 
/ = 1/2, for initial data such that = 0 and using a res¬ 
olution of Aa; = 0.003125 (see Appendix A for details on 
how to construct such initial data). For the initial data 
used here, according to lETHl a blow-up is expected at 
time T* = 9.98. The plots show both the initial data 
(dotted lines) and the numerical solution at time t = 10 
(solid lines), just after the expected blow-up. We clearly 
see how the eigenfield u> + has developed a large spike, 
while the lapse has developed a sharp gradient (the so¬ 
lution does not become infinite because the numerical 
method used has some inherent dissipation EH)- 

If one repeats the simulation at different resolutions 
one finds that the numerical solution converges up to a 
time t ~ 10, and after that convergence fails, indicating 
that even though the numerical solution continues past 
this time, we are no longer solving the original differential 
equations. To see this we can consider the convergence 
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Lapse function a Convergence factor 




FIG. 2: Evolution of the lapse function a for initial data rep¬ 
resenting a pulse traveling to the right. The dotted line shows 
the initial data, and the solid line the numerical solution at 
t = 10. Notice how a has developed a sharp gradient. 


of the constraint C a := D a — d x In a. Numerically this 
constraint will not vanish, but it should converge to zero 
as the resolution is increased. Define now the conver¬ 
gence factor as the ratio of the r.m.s norm of C a for a 
run at a given resolution and another run at twice the 
resolution. Figure 0 shows a plot of the convergence fac¬ 
tors as a function of time for runs done at five different 
resolutions: Ax = 0.05, 0.025, 0.0125, 0.00625, 0.003125. 
Since the numerical code is second order, the convergence 
factors should be close to 4. The figure shows that as the 
resolution is increased the convergence factors approach 
the expected value of 4 for t < 10 (the lines move up), 
but after that time they drop below 1 indicating loss of 
convergence. 

We have then found that blow-ups do develop for the 
Bona-Masso family of slicing conditions whenever (EDI 
does not hold, and we can even predict the precise time of 
blow up formation when f(a) is a constant. The question 
now is whether these blow ups are genuine shocks or not. 
To answer this question let us now consider the evolution 
of the eigenspeeds themselves along their corresponding 
characteristic lines. From we find that 

d t \± = ±d t (af 1 / 2 / g 1/2 ) , (3.13) 

and using now the evolution equations for a and g we 
obtain 


dt X±=a^(l-f-^pjK. (3.14) 

In a similar way we find for the spatial derivative 


— A^ 


D n 


«/' 


—jr ( / ^ n I n 


D n 


(3.15) 


FIG. 3: Convergence factors for the constraint C a for runs 
done at the five different resolutions: Ax = 0.05, 0.025, 
0.0125, 0.00625, 0.003125. As resolution increases, the con¬ 
vergence factors approach the expected value of 4 for t < 10 
(the curves move up), but after that time they drop sharply 
indicating loss of convergence. 


The last two equations together imply that 


dt\± + — 


, 1/2 


1 -/- 


a/' 


± s 1 ' 2 


D a 

T 


D n 


or in terms of the eigenfields 

dtX± + X±d x A^ = 

a/' 
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a 

1 J 2 


K =f 


D a 

/ 1/2 , 

(3.16) 




(3.17) 


If we now consider a region where = loq = 0, then 
the equation for A+ reduces to 


d t X + + A + 9a;A + — 0 . (3.18) 

This is nothing more than Burgers’ equation, the pro¬ 
totype for studying shock formation. It is easy to un¬ 
derstand how this equation implies genuine shocks: The 
equation shows that characteristic speeds are constant 
along their corresponding characteristic lines, so if these 
speeds where not uniform to begin with, and particularly 
if the derivative of A+ was initially negative at any point, 
the characteristic lines will inevitably cross. 

When the characteristics cross the spatial derivative of 
A + will become infinite, and as this derivative is given 
in terms of eigenfields, the eigenfields will blow up. This 
shows that the blow-ups we studied above correspond to 
places where the characteristic lines cross, i.e. they are 
genuine shocks. The use of the term “gauge shocks” to 
describe these pathologies is therefore entirely justified. 
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Position of characteristic lines 



Initial position 


Slices of Minkowski space-time 



FIG. 4: Positions of the characteristics during the simulation 
with respect to their initial positions, for a series of different 
times: t = 0, 2.5, 5, 7.5, 10. The characteristics were equally 
spaced a t = 0 (corresponding to a line at 45 degrees on the 
plot). By t = 10 this is no longer the case and a plateau has 
formed indicating that the characteristics are about to cross. 


FIG. 5: Initial hypersurface (dotted line) and final hypersur¬ 
face at t = 10 (solid line), as seen in Minkowski coordinates, 
for the simulation discussed in the text. The final slice has 
been moved back in time so that it lies on top of the initial 
slice at the boundaries. The initial slice is smooth, while the 
final slice has a sharp kink. 


Figure 0 shows the positions of the characteristics 
with respect to their initial positions for the simula¬ 
tion discussed above, at a series of different times: 
t = 0, 2.5, 5, 7.5,10. Notice that for the simulation con¬ 
sidered here / = 1/2 so the characteristic speed should 
be close to \f] ~ 0.707, we would then expect the lines on 
the plot to move up by approximately 0.707 x 10 = 7.07. 
From the figure we see that even though the characteris¬ 
tics were equally spaced a t = 0 (corresponding to a line 
at 45 degrees on the plot), at t = 10 this is no longer the 
case and a plateau has formed. This plateau indicates 
that those characteristics initially in the region around 
the origin are now all essentially at the same position 
x ~ 7, or in other words, they are about to cross. 

One could wonder what the gauge shocks imply about 
the geometry of the space-time being evolved. Since in 
this case we are simply evolving a foliation in Minkowski 
space-time it is clear that the background geometry re¬ 
mains perfectly regular, the only thing that can become 
pathological are the spatial hypersurfaces that determine 
the foliation. Figure O shows a comparison of the initial 
hypersurface and the final hypersurface at t = 10 as seen 
in Minkowski spacetime, using the data from the same 
numerical simulation discussed above (for an easier com¬ 
parison the final slice has been moved back in time so 
that it lies on top of the initial slice at the boundaries). 
The hypersurfaces are reconstructed by explicitly keep¬ 
ing track of the position of the normal observers in the 
original Minkowski coordinates during the evolution. No¬ 
tice how the initial slice is very smooth (it has in fact a 
Gaussian profile, see Appendix m, while the final slice 
has developed a sharp kink. This shows that the for¬ 


mation of a gauge shock indicates that the hypersurface, 
though still spacelike everywhere, is no longer smooth 
(its derivative is now discontinuous). 


IV. CONCLUSIONS 

I have studied the formation of gauge pathologies asso¬ 
ciated with the Bona-Masso family of slicing conditions 
in simple 1+1 relativity. The analysis shown here not 
only recovers previously known results about the exis¬ 
tence of such pathologies and the conditions necessary 
for avoiding them, but also proves that such pathologies 
arise from the crossing of the characteristic lines asso¬ 
ciated with the propagation of the gauge. This implies 
that such pathologies are genuine shocks and hence the 
name “gauge shocks” used to describe them is entirely 
justified. The gauge shocks appear as discontinuities in 
the lapse and spatial metric arising from smooth initial 
data, and represent places where the spatial hypersur¬ 
faces that determine the foliation of space-time become 
non-smooth. 

It is also important to point out that Ref. 0 introduces 
two different shock avoiding criteria called “indirect lin¬ 
ear degeneracy” and the “source criteria”. Neither crite¬ 
ria is fully understood at this point, and they are offered 
in Hj more as indicators that seem to work well in prac¬ 
tice than as rigorous conditions necessary for avoiding 
shocks. In this manuscript I decided against using either 
of those criteria precisely because they are not rigorous. 
The purpose here has been to show that for a simple case 
(1+1) one does not need to introduce any ad hoc criteria, 
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and one can instead show that gauge shocks develop by 
directly analyzing the evolution equations and showing 
that: 1) blow-ups in the eigenfields do happen (and can 
even be predicted in advance), and 2) the characteristic 
lines associated with the gauge propagation do cross. 


From this one can reconstruct the rescaled eigenfunction 
Q + and use Eq. (EH) to predict the time of shock for¬ 
mation given a specific form of h(x). 

In the simulation shown in Sec. m above, the profile 
function h was taken to be a simple Gaussian of the form 
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APPENDIX A 

In this Appendix I will discuss how to construct initial 
data that contains only waves propagating in one direc¬ 
tion. Notice that the simplest way to construct initial 
data would be to take a trivial initial slice in Minkowski 
coordinates, and use a non-trivial initial lapse function 
to introduce a later distortion. However, since in this 
case the initial extrinsic curvature vanishes, one sees 
from eh that inevitably modes traveling in both di¬ 
rection will be excited. We then conclude that in order 
to have waves propagating in only one direction one has 
to consider a non-trivial initial slice. 

The initial data considered here has already been dis¬ 
cussed in Q. and here I just present it for complete¬ 
ness. We start by considering an initial slice given in 
Minkowski coordinates (tjvf > %m') as 


t M = h(x M ) , (Al) 

with h a profile function that decays rapidly, like a Gaus¬ 
sian function. It is then not difficult to show that if we 
use x = Xm as our spatial coordinate, the spatial metric 
an extrinsic curvature turn out to be 


g = l-ti 2 =► D g = -2h’h"/g , (A2) 
K xx = —h"/sfg =» K = —h"/g. (A3) 

Assume that we want to have waves traveling to the 
right (the opposite situation can be done in an analogous 
way). This means that we want lo- = 0, which implies 

D a = y/f K = —y/f h"/g . (A4) 


This gives us the following differential equation a: 

d x a = —ay/f h"/ (l — h' 2 ) . (A5) 

In the particular case when / is a constant the above 
equation can be easily integrated to find 


1 - h '\ 77/2 
1 + h'J 


(A6) 


APPENDIX B 

Here I will briefly counter the arguments given by Bona 
et al. in @ to claim that gauge shocks are, in their words, 
completely discarded. Their main argument is based on 
looking at the so-called “foliation equation”, that looks at 
the foliation as the evolution of a scalar function T whose 
level sets correspond to the 3+1 slices. In jd] I showed 
that the Bona-Masso slicing condition can be written in 
covariant form as 

' g “" + ~ 75 *)) n “ n 

where is the unit normal vector to the hypersurfaces 

rig, = —aX7gT , (B2) 

and with the lapse function a given in terms of T as 

a = (—VgT V^T) -1/2 . (B3) 

Assume for a moment that / = 1. In that case Eq. m 
reduces to the standard wave equation. Consider now 
a specific point in spacetime and use locally flat coordi¬ 
nates. It is clear that in such coordinates the foliation 
equation has constant characteristic speeds equal to 1. 
This implies that the characteristic lines do not cross, 
and if the lines do not cross in some set of regular coor¬ 
dinates they won’t cross in any regular coordinates. We 
then conclude that in this case there are no shocks. 

But what happens when / ^ 1? Bona et al. argue that 
in that case one can always start with an arbitrary regular 
slice such that T = constant on that slice. From that 
one can easily see that, as long as the slice is spacelike, 
one can always recover all second partial derivatives of 
T from equation (IB III , which allows one to construct the 
next hypersurface. This is certainly true, but it does not 
follow from here that one can continue this procedure for 
any finite time, since all we have done is show that given 
regular initial data there is locally a regular solution of 
the foliation equation, but shocks are precisely solutions 
that fail after a finite time. 

In p] I showed that if one takes Eq. jnn> with / 7^ 1 
and considers locally flat coordinates, the characteristic 
speeds are not constant anymore. There is then no guar¬ 
antee that they will not cross after a finite time. In fact, 
in [H it is shown that if one applies the standard anal¬ 
ysis coming from the theory of PDE’s to Eq. (IB II) one 


V^Vj/T = 0 , (Bl) 
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finds that shocks will form unless condition (EH) holds. 
The problem with the argument of Bona et al. is there¬ 
fore that they failed to consider that when one deals with 
shocks, smooth initial data always guarantees smooth so¬ 
lutions locally , but not after a finite time. 

Bona et al. also present other arguments. For example, 
starting from their equation (6) for the speed of light 


dl 

c = — = ± 
dt 


dx l dxi 

7 « 1 THt =±a ’ <B4) 


they obtain the shock avoiding condition / = k/a 2 3 4 5 6 7 8 9 with 
k constant, which is clearly different from (13.71) . But 


in fact the last equality in the equation above is wrong, 
the speed of light is not equal to the lapse in a general 
coordinate system, one is missing a factor of \f r y™ (with 
x 1 the direction of propagation). Inserting this factor 
takes us back to the shock avoiding condition (1X71) . 


Bona et al. also consider cosmological-type scenarios 
(homogeneous and isotropic), and derive different con¬ 
ditions for avoiding blow-ups in the lapse (see equation 
(18) of f6]). But since in this case there is no propagation 
and just growth in place, the whole concept of shocks and 
characteristic speeds is simply non-applicable. 
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